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We report conductivity measurements of Cr-doped V 2 3 using a variable pres- 
sure technique. The critical behavior of the conductivity near the Mott-insulator 

o' 

to metal critical endpoint is investigated in detail as a function of pressure and 
O' 

temperature. The critical exponents are determined, as well as the scaling 

function associated with the equation of state. The universal properties of a 

liquid-gas transition are found. This is potentially a generic description of the 
O" 
$-* ; Mott critical endpoint in correlated electron materials. 
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Since the early recognition by Mott HJE1 that electron-electron interactions are responsible 
for the insulating character of many transition-metal oxides, extensive research over the last 
decade has demonstrated the key importance of this phenomenon for the physics of strongly 
correlated electron materials. Outstanding examples Q are superconducting cuprates, mangan- 
ites displaying colossal magnetoresistance, or fullerene compounds. There are two routes for 
achieving a metallic state, starting from a Mott insulating material. The first is to introduce 
charge carriers by doping. The second, closely connected to Mott's original ideas, is to reduce 
the ratio U/W between the typical strength of local Coulomb repulsion (U) and the typical 
kinetic energy of the relevant electrons (W). This can be achieved in practice, in some materi- 
als, by selected atomic substitutions or by applying pressure. The most widely studied example 
ElEl SHOD is Cr-doped Vanadium sesquioxide (V\- x Cr z ) 2 03 which displays a transition from 
a paramagnetic Mott insulator to a strongly correlated metal. The transition into the metallic 
state can be triggered by lowering temperature (at sufficiently small chromium concentration 
x), by decreasing x or by increasing pressure P (early studieslSJ |3 have revealed that 
decreasing concentration by Aa; ~ —0.01 corresponds to an applied 

pressure of AP ~ 4kbar). The transition is first-order, with a significant reduction of the lat- 
tice spacing through the insulator-to-metal transition, indicating a coupling between electronic 
and lattice degrees of freedom. The first-order transition line in the (P, T) - or (x, T) -phase 
diagram ends in a second-order critical endpoint (P c , T c ). 

We report on transport experiments which allow for a precise identification of the critical 
behavior associated with this critical endpoint, a question of fundamental importance in under- 
standing the Mott transition. Recent theoretical developments have proposed a description of 
the critical behavior in simplified purely electronic models (and also of the crossovers between 
distinct transport regimes close to the critical point). Despite extensive studies of this material, 
the critical behavior has not been elucidated so far experimentally. The key technique used in 



the present work is to perform conductivity measurements as a function of continuously vary- 
ing pressure, at constant temperature (see ref. Q for an early study). This is particularly well 
adapted to the present situation in which the transition line in the (T, P) plane is very sharp. In 
contrast, the work of Kuwamoto et al. JE) investigated conductivity as a function of temperature 
for a discrete set of chromium concentrations x. 

We use an isopentane liquid pressure cell, and control the value of pressure with an accuracy 
of 1 bar. Conductivity is measured at constant regulated temperature with an accuracy of order 
0.1 K, as a function of pressure, using a standard four-probe method. All our measurements 
were performed on crystals of (V .989 Cr .on)20 3 grown using the skull-melter technique fol- 
lowed by appropriate annealing [9|. The choice of a Cr concentration x = 0.011 is such that 
the sample is on the insulating side of the transition at ambient pressure but that a moderate 
pressure (a few kilobars) drives the system into the metallic state (or, alternatively, a decrease 
in temperature). This is visible on the data set (Fig. 1A), which displays the conductivity cr as a 
function of pressure P, for several temperatures in the range 290K < T < 485K. These data are 
obtained by decreasing pressure from P = 6 kbar down to ambient pressure, going from a high- 
conductivity metallic regime to a low-conductivity insulating regime. For temperatures smaller 
than the critical temperature T c this transition is discontinuous, with a sudden jump of the con- 
ductivity. In order to locate precisely this critical point and to demonstrate the first-order nature 
of the transition, we have performed hysteresis experiments in which the conductivity is mea- 
sured during increasing and decreasing pressure sweeps at a slow rate of order 25 bar/min (Fig. 
1A). From the difference between the measured conductivities in these two sweeps (Fig. S2), 
two characteristic pressures can be identified, Pm(T) and Pi(T) (P M < Pj), corresponding 
respectively to the lowest pressure at which a metallic state can be sustained while decreasing 
pressure (Pa/), and to the highest pressure at which an insulating state can be sustained while 
increasing pressure (Pi). These two spinodal lines , plotted as a function of temperature on Fig. 



IB, merge at the critical endpoint (P c , T c ). We can then estimate: P c ~ 3738 bar, T c ~ 457.5 K. 
Varying pressure rather than temperature is essential for a precise determination of T c , which 
is compatible with the early estimate (~ 450K) by Kuwamoto et al.[8|. At the critical tem- 
perature, the pressure dependence of a(P, T c ) becomes singular, with a vertical tangent at the 
critical pressure P = P c (Fig. 2A). For T > T c , this singular behavior is replaced by a contin- 
uous variation of the conductivity with pressure, which nevertheless defines a sharp crossover 
line in the (P, T) phase diagram (as also depicted in Fig. IB). This crossover line extrapolates 
to a temperature of order ~ 500K for the pressure (~ 5kbar) corresponding to the pure V2O3 
compound. Interestingly, the location of this crossover coincides with the one detected in early 
NMR experiments! 10| . 

We now show that the critical singularities found in the vicinity of the critical endpoint 
(P c , T c ) can be analyzed in the framework of the scaling theory of the liquid-gas transition of 
classical systems[l 1J. The analogy between the latter and the finite-temperature Mott transition 
has been emphasized early on by Castellani et al.WlA (see also Ref. [5|). The insulating phase 
(in which the Vanadium is mainly in the V 3+ state, corresponding to the d 2 configuration) can 
be pictured as a "gas" phase with a low density of double occupancies or holes (corresponding 
to V 2+ and V 4+ , or d 3 and d 1 , respectively). The metallic phase corresponds to a "liquid" 
with a sizeable density of holes and double occupancies. Recently, this analogy has been given 
firm theoretical foundations within the framework of a Landau theory lfT3l fT4ll derived from 
dynamical mean-field theory (DMFT) lfT5l . In this framework, a scalar order parameter <\> is 
associated with the low-energy electronic degrees of freedom which build up the quasiparticle 
resonance in the strongly correlated metallic phase close to the transition. This order parameter 
couples to the singular part of the double occupancy (hence providing a connection to the picture 
described above), as well as to other observables such as the Drude weight or dc-conductivity. 
Because of the scalar nature of the order parameter, the transition falls in the Ising universality 



class. Coupling to lattice degrees of freedom can also be included in the theory iflol without 
changing this conclusion. In the following we denote by r the scaling variable corresponding 
to the temperature scaling axis in the Ising model analysis (i.e. to the term rcj) 2 in the Landau 
functional) and by h the scaling axis corresponding to magnetic field (i.e. to the symmetry- 
breaking term —h(p). These scaling variables are a priori linear combinations of (T — T c )/T c 
and (P — P c )/P c . However, our data are compatible with no or little mixing, so that we choose 

in all the following: r = (T - T c )/T c -\ , h = (P - P c )/P c -\ (the dots indicate higher 

order terms). Denoting by a c = a(P c ,T c ) (~ 449.5 Cl~ 1 cm~ 1 ) the measured conductivity at 
the critical point, it is expected that a(P, T) — a c depends linearly on the order parameter (0) 
close to the critical point. (This can be explicitly proven in the context of DMFT). At T = T c , 
this implies a critical singularity of the form: a(P, T c ) — a c ~ h 1 ^ 5 with 5 the critical exponent 
associated with the singular dependence of the magnetization at the critical point in the Ising 
model. The data in Fig. 2A are very well fitted by this form, as demonstrated in the inset. 
Over more than two decades in h oc (P — P c )/P c , we find the best-fit value of the exponent 
to be 5 ~ 3, i.e. the mean-field value. In a narrow pressure interval (AP ~ 10 bar) close to 
the critical pressure, indication for a crossover towards a value 5 ~ 5 is found, close to the 
three-dimensional (3D) Ising value 5 ~ 4.814. 

We now address the critical behavior away from T c by studying the temperature dependence 
of the conductivity in the following manner. For T < T c , we focus on the conductivity of the 
metallic state, at the high-pressure boundary of the coexistence region. That is, we consider 
a*(T) = a met (Pi(T),T) — a c with Pi(T) the spinodal of the insulating phase. This quantity, 
plotted in Fig. 2B, is expected to display the critical behavior of the order parameter, by analogy 
with the liquid-gas transition: a*(T) ~ (— r)' 3 with r oc (T — T c )/T c . As shown in the inset 
(Fig. 2B), a mean-field value of the critical exponent /? ~ 0.5 is found to fit the data over 
almost two decades away from the critical point. In a narrow temperature interval AT ~ 4K 



close to T c (AT/T C ~ 0.01), some indication for a crossover towards a non mean-field value 
(3 ~ 0.34 is found, close to the 3D Ising value (3 ~ 0.327. We also consider the derivative 
of the conductivity with respect to pressure, in the metallic state, taken on the same spinodal 
line: x(T) = (da met (P,T)/dP)\ P=Pl i yT y This quantity can be defined as well for T > T c 
by taking the derivative at the inflection point of the cr(P) curve (Fig. 1A). Following the 
liquid-gas analogy, it corresponds to the magnetic susceptibility in the equivalent Ising model: 
X oc d(4>)/dh. As shown in Fig. 2C, it is found to diverge as x ~ X+/(T — T c ) y for T > T c and 
as x ~ X-/(T C — T) 7 ' for T < T c . The exponent 7, as well as the (universal) amplitude ratio 
X+/x—> are found to be close to their mean-field values: 7 = 7' = !, x+/X- — 2. Very close to 
T c , the noise in the numerical derivative involved in the determination of x prevents a reliable 
determination of deviations from mean-field, in contrast to the above study of the conductivity 
itself. 

Finally, we demonstrate that the whole set of conductivity data in the metallic phase can be 
scaled onto a universal form, which can be written as: 

a met (P,T) -a c = {5h) l/S U (|^=iy) (D 

In this expression, 5h = h — hi denotes the difference between the "field" h = (P — P c )/P c 
and its value on the spinodal line of the insulating phase hj = (Pi — P c )/P c , i.e.: Sh = 
(P — Pi(T))/P c . This amounts to a simple shift of the field variable on the standard form ifTTIl 
of the universal equation of states near a liquid-gas critical point. The functions / + and /_ are 
universal scaling functions which apply for T > T c (r > 0) and T < T c (r < 0), respectively. 
When written in this form, the equation of state is such that the order parameter <r*(T) defined 
above is recovered when the limit Sh — » is taken in the right-hand side of Eq. 1 . The pressure- 
dependent data sets for many different temperatures have been plotted in this manner (Fig. 3), 
in which the two exponents 7 and 5 were taken as adjustable parameters (Fig. S3) in order to 



obtain the best collapse of all the data points onto single curves. This leads to values of these 
exponents close to the mean-field ones 7 ~ 1, 5 ~ 3, which provides a strong check on the 
individual determination of each critical exponent performed above. The quality of the scaling 
is seen to be excellent over a very large range of variation of the scaling variables (several 
decades). It is apparent that the scaling functions obey the expected asymptotic behaviors: 
f+(x <C 1) ~ x 1-1 / 5 , f-(x <1)~ x^ 1 / 6 and f±(x ^> 1) ~ const.. This finding is essential to 
ensure that Eq. 1 be compatible with the critical behavior of the order parameter a* at small and 
large field, for both T < T c and T > T c , investigated previously in Fig. 2. It also implies that 
the critical exponents obey the relation 7 = (3(5 — 1), in agreement with the above determination 
of 13. 

These universal scaling properties of the pressure- and temperature-dependent conductivity 
experimentally demonstrate that the electronic degrees of freedom undergo a liquid-gas phase 
transition at the Mott critical endpoint. Critical exponents and a universal scaling function have 
been determined. Our results are consistent with mean-field values over most of the parameter 
range, with some indication for three-dimensional Ising behavior very close to the transition. A 
possible explanation for why the range of validity of mean-field theory is so large can be put 
forward by analogy with the theory of conventional superconductors. There, the key point is the 
existence of a very large length scale (the pair coherence length), much larger than the lattice 
spacing (or the Fermi wavelength). Here, the Mott insulator can be thought of as a state in which 
holes and doubly occupied sites form bound states due to their Coulomb interaction. The spatial 
extension £ of these bound states is related to their energy (the Mott gap A) by A ~ h 2 /(2m^ 2 ). 
Given the measured value of A in samples close to the transition, this leads to the conclusion 
that £ is indeed a large length-scale, of order a few nanometers. Finally, we emphasize that our 
results provide experimental support to the early idea of Ref. lfT2l and to recent theories of the 
Mott critical endpoint based on the dynamical mean-field (DMFT) approach! 13l fl4lfT5l . While 



further effort should be devoted to the inclusion of lattice degrees of freedom in these theories, 
simplified treatments of these effects ifTol do emphasize the key role of electronic degrees of 
freedom in the transition. 
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Figure captions 

Fig. 1A: Conductivity as a function of decreasing pressure, for temperatures ranging from 
T = 485K (> T c = 457.5K, orange curves) down to T = 290K (< T c , blue curves). The dark 
yellow curve is the conductivity at T c . Only a selected set of values of T has been displayed, for 
clarity. (For a two-dimensional plot of the data, see Fig. SI) Examples of an hysteresis cycle 
are shown for T = 290K and T = 348K. For a plot of the difference of conductivities measured 
in decreasing and increasing pressure sweeps, see Fig. S2. 

Fig. IB: Phase diagram of Cr-doped Vanadium Sesquioxide Vx_ x . Cr^Os as a function of 
pressure and temperature, in the range lbar < P < 6kbar and 290K < T < 500K investigated 
in this work. At a given temperature T, the metallic state can be obtained for pressures higher 
than the spinodal pressure Pm(T), and the insulating state for pressures lower than the spinodal 
pressure Pi(T). These two spinodal lines delimit a pressure range Pm < P < Pi in which the 
two states coexist (hatched region on the figure). This coexistence region closes at the critical 
endpoint (P C ,T C ) (P c ~ 3738bar, T c ~ 457.5K). The crossover line above this point (dashed) 
corresponds to the inflection point in the cr(P) curves. 

Fig. 1C: Schematic global phase diagram of Cr-doped Vanadium Sesquioxide Vi- X Cr^Os as 
a function of pressure and temperature, deduced from Ref. Q. 

Fig. 2 These plots demonstrate how the critical exponents 5, (3 and 7 can be inferred from the 

study of the conductivity and of its derivative with respect to pressure (see text). 

Fig. 2A: At the critical temperature T = T c , the conductivity a is plotted as a function of 

pressure. The (plain) red line is a fit to a — a c ~ (P — P c ) 1 ^ ', with 6 = 3. The use of a 

logarithmic scale (inset) confirms this value, and also reveals a non mean-field regime for P 

close to P c . 

Fig. 2B: Order parameter a*(T) = a(P I (T),T) - o c vs. T/T c , for T < T c . The line is a fit to 
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(T c — T) 13 with (3 = 0.5. The inset (logarithmic scale) reveals a non mean-field regime close to 
T c . 

Fig. 2C: Derivative of the conductivity (analogous to a susceptibility %, as described in text), 
for T < T c and T > T c . The plain lines are fits to x±\T - T c \~i , with 7 = 1 and x+/X- = 2. 
Fig. 3: Scaling plot of the conductivity onto a universal equation of state. The whole data set 
in the metallic state has been used in order to plot (a — a c )/(P — Pi) 1 ^ 5 vs. (P — Pj)/(T — 
r p^ps/{8-\) ^ as d escr ibed in the text. The data collapse onto two universal curves for T > T c and 
T < T c , corresponding to the universal scaling functions f± in Eq.(l). 
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